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ABSTRACT 

The inability of standard models to explain the flux ratios in many 4-image gravitational lens systems 
has been presented as evidence for significant small-scale structure in lens galaxies. That claim has 
generally relied on detailed lens modeling, so it is both model dependent and somewhat difficult 
to interpret. We present a more robust and generic method for identifying lenses with small-scale 
structure. For a close triplet of images created when the source lies near an ideal cusp catastrophe, 
the sum of the signed magnifications should exactly vanish, independent of any global properties of 
the lens potential. For realistic cusps, the magnification sum vanishes only approximately, but we 
show that it is possible to place strong upper bounds on the degree to which the magnification sum 
can deviate from zero. Lenses with flux ratio "anomalies," or fluxes that significantly violate the 
upper bounds, can be said with high confidence to have structure in the lens potential on scales of the 
image separation or smaller. Five observed lenses have such flux ratio anomalies: B2045+265 has a 
strong anomaly at both radio and optical/near-IR wavelengths; B0712+472 has a strong anomaly at 
optical/near- IR wavelengths and a marginal anomaly at radio wavelengths; 1RXS J1131— 1231 has a 
strong anomaly at optical wavelengths; RX J0911+0551 appears to have an anomaly at optical/near- 
IR wavelengths, although the conclusion in this particular lens is subject to uncertainties in the typical 
strength of octopole density perturbations in early-type galaxies; and finally, SDSS J0924+0219 has 
a strong anomaly at optical wavelengths. Interestingly, analysis of the cusp relation does not reveal 
a significant anomaly in B1422+231, even though this lens is known to be anomalous from detailed 
modeling. Methods that are more sophisticated (and less generic) than the cusp relation may therefore 
be necessary to uncover flux ratio anomalies in some systems. Although these flux ratio anomalies 
might represent either milli-lensing or micro-lensing, we cannot identify the cause of the anomalies 
using only broad-band flux ratios in individual lenses. Rather, the conclusion we can draw is that 
the lenses have significant structure in the lens potential on scales comparable to or smaller than the 
separation between the images. Additional arguments must be invoked to specify the nature of this 
small-scale structure. 

Subject headings: cosmology: theory — dark matter — galaxies: formation — gravitational lensing 
— large-scale structure of universe 



1. INTRODUCTION 

Gravitational lens modeling has had remarkable suc- 
cess handling increasingly precise measurements (e.g., 
Barkana et al. 1999; Patnaik et al. 1999; Trotter, Winn, 
& Hewitt 2000) and increasingly sophisticated datasets 
including Einstein ring images (Keeton et al. 2000; 
Kochanek, Keeton, & McLeod 2001) and/or stellar dy- 
namical data (Romanowsky & Kochanek 1999; Koop- 
mans & Treu 2002, 2003; Treu & Koopmans 2002, 2003; 
Koopmans et al. 2003b). Lens modeling has even clar- 
ified the properties of complex systems with more than 
one lens galaxy and/or more than one background source 
(Cohn et al. 2001; Rusin et al. 2001; Keeton & Winn 
2003; Koopmans et al. 2003b). However, the notable 
success has largely been restricted to the number and 
configuration of lensed images. The flux ratios between 
the images, at least in lenses with four or more images, 6 
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have long resisted explanation. 

Until recently the persistent problem with flux ratios in 
4-image lenses (e.g., Kent & Falco 1988; Falco, Lehar, & 
Shapiro 1997; Keeton, Kochanek, & Seljak 1997) received 
little attention, perhaps because the number of 4-image 
lenses was relatively small, and because it seemed pos- 
sible to appeal to electromagnetic — non-gravitational 
- effects such as extinction by dust or scattering by 
hot gas. However, the number of lenses with apparently 
anomalous flux ratios is growing rapidly (e.g., Inada et 
al. 2003; Sluse et al. 2003; Wisotzki et al. 2003). More- 
over, direct evidence suggests that electromagnetic ef- 
fects, while present in some lenses, cannot explain most 
of the anomalies (Falco et al. 1999; Winn et al. 2001, 
2002; Koopmans et al. 2003a). The problem with flux 
ratios therefore appears to be real. 

It also turns out to have interesting and important im- 
plications for astrophysics and cosmology. When Mao 
& Schneider (1998) made the first systematic analysis of 
the flux ratio problem, they realized that the anomalies 
might be attributed to gravitational effects omitted from 
standard lens models, namely small-scale structure in the 
lens galaxy. The key insight was that since flux ratios 
are determined by second derivatives of the lens poten- 

models. 
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tial, they are much more sensitive to small-scale structure 
than the image positions (which are determined by first 
derivatives of the potential); so models that lack small- 
scale structure might successfully reproduce the image 
positions but fail to fit the flux ratios. 

One possible source of small-scale structure is clumps 
of dark matter of mass ~ 10 6 -10 9 M Q left over from the 
hierarchical galaxy formation process in the Cold Dark 
Matter (CDM) paradigm (Metcalf & Madau 2001; Chiba 
2002; Dalai & Kochanek 2002). This possibility has gen- 
erated significant interest because it relates to current 
questions about the validity of CDM on small scales. 
The discrepancy between the predicted abundance of 
dark matter clumps and the observed abundance of dwarf 
galaxy satellites around the Milky Way has been inter- 
preted as a fundamental problem with CDM (Klypin et 
al. 1999; Moore et al. 1999), which may signal a need for 
new physics for the dark matter (e.g., Spergel & Stein- 
hardt 2000; Colin, Avila-Reese, & Valenzuela 2000; Hu, 
Barkana, & Gruzinov 2000). Alternately, the discrep- 
ancy may simply indicate poor understanding of the as- 
trophysical processes that determine whether or not a 
clump of dark matter hosts a visible dwarf galaxy (Bul- 
lock, Kravtsov, & Weinberg 2000; Benson et al. 2002; 
Somerville 2002; Stoehr et al. 2002; Hayashi et al. 2002). 
If lens flux ratios can be used to probe dark matter 
clumps, that will provide the cleanest way to distinguish 
these two very different hypotheses, and more generally 
to resolve the controversy about whether CDM does or 
does not over-predict small-scale structure (e.g., Flores & 
Primack 1994; Moore 1994; Spergel & Steinhardt 2000; 
Debattista & Sellwood 2000; de Blok, McGaugh, & Ru- 
bin 2001; Keeton 2001a; van den Bosch & Swaters 2001; 
Weiner, Sellwood, & Williams 2001; de Blok & Bosma 
2002; Kochanek 2003). Early results indicate that the 
statistics of flux ratio anomalies imply a clump popula- 
tion that agrees well with CDM predictions and validates 
cold dark matter (Dalai & Kochanek 2002; Kochanek & 
Dalai 2003), but the importance of the conclusion de- 
mands further study. 

A second interesting possibility is that the small-scale 
structure implied by flux ratio anomalies is simply stars 
in the lens galaxy (Chang & Refsdal 1979; Irwin et al. 
1989; Wozniak et al. 2000; Schechter & Wambsganss 
2002). In this case, flux ratio anomalies offer a unique 
probe of the relative contributions of stars and dark mat- 
ter to the surface mass density at the image positions 
(Schechter & Wambsganss 2002), which would be inter- 
esting because the amount of dark matter contained in 
the inner regions of elliptical galaxies is still not well 
known (e.g., Gerhard et al. 2001; Keeton 2001a; Bor- 
riello, Salucci, & Danese 2003; Rusin, Kochanek, & Kee- 
ton 2003b). Yet a third possibility is that the small- 
scale structure is not localized like dark matter clumps 
or stars, but is more global like small disk components in 
bulge-dominated systems, Fourier mode density fluctu- 
ations, tidal streams, etc. (e.g., Mao & Schneider 1998; 
Evans & Witt 2002; Quadri, Moller, & Natarajan 2003; 
Moller, Hewett, & Blain 2003). If this is the case, then 
lensing can be used to search for such structures whether 
they are traced by the luminous components of galaxies 
or not. 

These three disparate applications all rest on a com- 
mon foundation: the identification of lenses with flux 



ratio anomalies that indicate small-scale structure. That 
identification is most unambiguous when time variability 
(e.g., Wozniak et al. 2000; Schechter et al. 2003) or re- 
solved spectra of the images (e.g., Moustakas & Metcalf 
2003; Wisotzki et al. 2003) clearly indicate microlensing 
by stars in the lens galaxy, or when the resolved shapes 
of the images indicate structure on the scale of dark mat- 
ter clumps (e.g., Metcalf 2002). Until such data become 
available for the majority of lenses, however, we need 
a method to identify anomalies using only broad-band 
flux ratios. Besides, such a method will be needed to se- 
lect candidates for the expensive follow-up observations 
(monitoring, spectroscopy, or high-resolution imaging). 

To date, the usual approach has been to use detailed 
lens modeling to interpret broad-band flux ratios and 
draw conclusions about, for example, the abundance of 
dark matter clumps (e.g., Dalai & Kochanek 2002; Met- 
calf k, Zhao 2002; Kochanek & Dalai 2003). This ap- 
proach is vulnerable to the criticism that the results de- 
pend on the sorts of lens potentials used in the modeling. 
The argument has two parts. First, many of the com- 
monly used families of lens potentials implicitly possess 
global symmetries, which lead to invariant magnification 
relations that are "global" in the sense that they involve 
all four images (Dalai 1998; Witt & Mao 2000; Dalai 
& Rabin 2001; Hunter & Evans 2001; Evans & Hunter 
2002). If a fit is poor because the data fail to satisfy these 
relations, that does not automatically constitute a flux 
ratio anomaly; it may simply indicate that the assumed 
relations are too restrictive, and that small, unremark- 
able deviations from the assumed symmetries are needed. 
The conceptual difficulty here is that one is trying to use 
global relations to draw conclusions about structure on 
smaller, more local scales. The second part of the argu- 
ment is that there is a large difference in scale between 
the image separations (~0'.'2-2") and the scales relevant 
for dark matter clumps (~10~ 3 ") or stars (~10 -6 "). If 
the flux ratio anomalies are in fact due to structures that 
are intermediate between these scales, then they may not 
necessarily imply the presence of dark matter clumps or 
stars (Evans & Witt 2002; Quadri et al. 2003; Moller et 
al. 2003). 

To address the first part of the criticism, we seek a 
method of identifying flux ratio anomalies that is local 
rather than global, i.e., a method that is sensitive only 
to structures smaller than the scales probed by the im- 
age positions. Fortunately, one can do this by appealing 
to simple, generic relations between the image magnifi- 
cations that should be satisfied for images in "fold" or 
"cusp" configurations (defined in The magnification 
relations are derived from local properties of the lens 
mapping and are in principle independent of the global 
mass model. They can be violated only if there is signifi- 
cant structure in the lens potential in scales smaller than 
the separations between the images (see Mao & Schnei- 
der 1998). In practice, however, the situation is compli- 
cated by the fact that the caustics in real lens systems 
only approximate ideal folds and cusps in some low-order 
expansion of the potential near the critical point; higher- 
order terms introduce deviations from the fold and cusp 
geometries. Real lenses therefore need to obey the ideal 
magnification relations only approximately. Because the 
accuracy with which the relations should hold depends 
on the distance of the images from the critical point and 
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on properties of the lens potential, it is not straightfor- 
ward to judge a priori the significance of an apparent 
violation. 

Our goal is to understand the magnification relations 
in realistic lens potentials and to determine how well they 
can be used to identify flux ratio anomalies. In this paper 
we focus on cusp configurations, because as the highest 
order stable singularities in lensing maps (see Schneider, 
Ehlcrs, & Falco 1992; Petters, Levine, & Wambsganss 
2001) cusps are amenable to analytic study, and cusp 
configurations are easy to identify. 7 We will address fold 
configurations in subsequent work. We study the degree 
to which the ideal cusp relation can be violated due to 
various properties of the lens potential: the radial den- 
sity profile, ellipticity, and multipole density perturba- 
tions of the lens galaxy, and the external tidal shear from 
the lens environment. Using both analytic and numer- 
ical methods we derive upper bounds on the deviation 
from the ideal cusp relation for realistic lens potentials 
that lack significant small-scale structure. We then argue 
that finding larger deviations in observed lenses robustly 
reveals flux ratio anomalies and indicates the presence of 
some sort of small-scale structure. 

We assert that, even though we adopt specific fami- 
lies of lens potentials, our analysis is more general than 
explicit modeling. One reason is that we have a bet- 
ter distinction between global and local properties of the 
lens potential. For example, a global m = 1 mode (i.e., 
non-reflection symmetry) would affect conclusions about 
anomalies in direct modeling, but not in our analysis. 
A second reason is that we consider quite general forms 
for the lens potential and take care to understand which 
generic features affect the cusp relation. A third point is 
that our results are less modeling dependent, less subject 
to the intricacies of fitting data and using minimization 
routines. A fourth advantage of our analysis is that, 
rather than simply showing that standard models fail 
to fit a lens, it clearly diagnoses why. We believe that 
these benefits go a long way toward establishing that 
small-scale structure in lens galaxies is real and can be 
understood. 

We must address a question that is purely semantic but 
nevertheless important: Where do we draw the line be- 
tween a normal "smooth" lens potential and "small-scale 
structure"? Taking a pragmatic approach, we consider 
"smooth" to mean any features known to be common 
in (early-type) galaxies: certain radial density profiles, 
reasonable ellipticities, small octopole modes represent- 
ing "disky" or "boxy" isophotes, and reasonable exter- 
nal shears. We consider "small-scale structure" to be 
anything whose presence in early-type galaxies would be 
notable. This can include stars — although stars are ob- 
viously abundant in galaxies, detecting the gravitational 
effects of individual stars is still interesting — and dark 
matter clumps, which seem to have generated the most 
interest. But it may also include tidal streams, massive 
or offset disk components (see Quadri et al. 2003; Moller 
et al. 2003), large-amplitude multipole density fluctua- 
tions (see Evans & Witt 2003), etc. We emphasize that 
our analysis, or indeed any analysis that considers only 

7 A close triplet of images always indicates a cusp configuration; 
but a close pair of images could be associated with either a fold or 
a cusp. 



the image positions and broad-band flux ratios in in- 
dividual lenses, cannot distinguish between these types 
of small-scale structure. The most general conclusion 
we can draw from flux ratio anomalies is that the lens 
potential contains structure on scales comparable to or 
smaller than the separation between the images. Further 
data and analysis is required to determine the nature of 
the small-scale structure (e.g., Wozniak et al. 2000; Met- 
calf 2002; Kochanek & Dalai 2003; Moustakas & Metcalf 
2003; Schechter et al. 2003; Wisotzki et al. 2003). 

The layout of the paper is as follows. We begin in 
by reviewing quadruple imaging and introducing a way 
to characterize 4-image configurations quantitatively. (In 
this paper we consider only 4-image lenses.) In fJ3]we dis- 
cuss cusp image configurations and present the generic, 
universal relation that should be obeyed by the image 
magnifications for sources near an ideal cusp. We then 
test this ideal relation, first using analytic results for sim- 
ple lens potentials (SQJ, and then with Monte Carlo sim- 
ulations of realistic lens populations (jJHJl. In jJJ]] we ap- 
ply the cusp relation to observed lenses, using violations 
of the relation to identify lenses that require small-scale 
structure. We offer our conclusions in jJJ| Several ap- 
pendices present supporting technical material. In Ap- 
pendix [S] we derive the universal relations between the 
image positions and magnifications for sources near an 
ideal cusp. In Appendix 151 we obtain exact analytic so- 
lutions of the lens equation for two families of lens po- 
tentials, which can be used to obtain exact analytic ex- 
pressions for the realistic cusp relation. 

2. CHARACTERIZING 4-IMAGE LENSES 

Nineteen quadruply-imaged lens systems have ap- 
peared in the literature, and they are listed in Table ITTT1 
This count includes only systems that have exactly four 
images of a given source, and where the images appear 
point-like at some wavelength. It includes the 10-image 
system B1933+503, which is complex only because there 
are three distinct sources; none of the sources has an im- 
age multiplicity larger than four (Sykes et al. 1998). By 
contrast, it excludes PMN J0134-0931 and B1359+154 
because they have multiplicities larger than four due to 
the presence of multiple lens galaxies (Rusin et al. 2001; 
Keeton & Winn 2003; Winn et al. 2003). One other lens, 
0047—2808, is almost certainly quadruply-imaged as well 
(Warren et al. 1996, 1999; Koopmans & Treu 2003), but 
its lack of point-like images makes it difficult to analyze 
with the usual techniques used for point-like systems. 

Mathematically, quadruple imaging can be described 
in terms of the critical curves and caustics of the lens 
potential. (See the monographs by Schneider et al. 1992 
and Petters et al. 2001 for thorough reviews of lens the- 
ory.) Critical curves are curves in the image plane where 
the lensing magnification is formally infinite, and caustics 
are the corresponding curves in the light source plane. 
The properties of these curves can be studied with catas- 
trophe theory; for our purpose the important result is 
that the astroid-shaped caustic that is associated with 
quadruple imaging has a generic shape that leads to three 
generic configurations of 4-image lenses (see Figure ^| . 
Sources near a cusp in the caustic produce "cusp" con- 
figurations with three of the images lying close together 
on one side of the lens galaxy. Source near the caustic 
but not near a cusp produce "fold" configurations with 
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Fig. 1. — The three basic configurations of 4-image lenses: 
fold (top), cusp (middle), and cross (bottom). In each panel, the 
figure on the left shows the caustics and source position in the light 
source plane, while the figure on the right show the critical curves 
and image positions in the image plane. 
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two of the images lying close together. Sources not close 
to the caustic produce relatively symmetric "cross" con- 
figurations. 

Although it may seem easy to label an observed lens 
as a fold, cusp, or cross, the categories actually blend 
together so it is important to develop a more quantitative 
way to characterize image configurations. To quantify a 
triplet of images (as in a cusp configuration), let d be 
the maximum separation between the three images, and 
let 9 be the opening angle of the polygon spanned by 
the three images, measured from the position of the lens 
galaxy. Each 4-image lens has four distinct triplets and 
hence four values of 9 and d. We can identify image 
triplets associated with cusps as those where 9 and/or d 
is small (see Figure^. Even though there is no rigorous 
definition of when 9 and d are "small" enough to indicate 
a cusp configuration, we shall see below that these are 
useful quantities for characterizing the range of image 
configurations. 

3. UNIVERSAL MAGNIFICATION RELATION FOR CUSPS 

In this section we briefly review the lensing of a source 
close to and inside an ideal cusp and present the magni- 
fication relation used in our analysis. Appendix 1X1 dis- 
cusses lensing near a cusp in considerably more detail, 
and presents additional position and magnification rela- 
tions for cusp images. This analysis applies to ordinary 
cusps; it may not be valid for so-called ramphoid cusps 
(or cusps of the second kind), but such cusps have not 
been observed and are expected to be rare in lensing sit- 
uations of astrophysical interest (see Petters & Wicklin 
1995; Oguri et al. 2003). 

In the vicinity of a cusp, the lens equation relating the 
source position u to the image position 9 can be written 



Fig. 2. — Sample image triplets for the image configurations 
from Figure^ together with the values of the opening angle 8 and 
the separation d (in units of the Einstein radius of the lens). The 
star shows the position of the lens galaxy in each system. 



to third order in 8 as a polynomial mapping, 



'2 ' 



u 2 = b 9i 9 2 



(1) 



The coordinates u and 6 are local orthogonal coordinates 
that are related to the global coordinates y and x of the 
lens system by u = My and 6 = Mx, where the transfor- 
mation matrix M depends on the lens potential. For the 
simple cases that we study in S0]and Appendix iBl M is 
the identity matrix and the 6 and u coordinate systems 
are simply the x and y coordinate systems translated so 
the cusp point is at the origin. The constant coefficients 
a, &, and c are given by derivatives of the potential at 
the critical point (see eq. IA6I in Appendix [X I . 

Solving for 9\ in the left-hand side of eq. and sub- 
stituting into the right-hand side, one obtains a cubic 
equation for 9i that depends on a, 6, c, and the source 
position u. Inside the caustic, there are three real so- 
lutions to this cubic equation, and thus three images of 
the source. It is possible to derive six independent rela- 
tions between the positions and magnifications of these 
images. Unfortunately, only one of these relations can 
be recast to depend only on directly observable proper- 
ties: the well-known magnification sum rule (Schneider & 
Weiss 1992; Zakharov 1995; Petters et al. 2001, p. 339), 



Ml + A*2 + A*3 = , 



(2) 



where the Hi are the signed magnifications of the three 
images. The other relations depend on properties that 
are not directly observable, such as the position of the 
source or the mapping coefficients a, &, and c. 
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4. THE CUSP RELATION IN SIMPLE LENS POTENTIALS 

The derivation of the ideal cusp relation eq. @ relics 
on the assumption that the lensing map has the polyno- 
mial form of eq. • Since this form is a truncated Taylor 
series expansion near the cusp point, we should expect 
the cusp relation to be exact only for sources asymp- 
totically close to the cusp. In this section we begin to 
quantify the deviation from the ideal cusp relation that 
arise from the higher order terms in the lensing map, us- 
ing simple examples to illustrate the effects of the radial 
profile, cllipticity, shear, and multipole perturbations of 
the lens potential. 

The magnifications appearing in the cusp relation are 
not directly observable, but we can follow Mao & Schnei- 
der (1998) and divide out the unknown source flux by 
defining the dimensionless quantity 

d +M2 + M3I = \F!+F 2 + F 3 \ 

CUSP " | Ml | + | M2 |+| M3 | |iq| + |F 2 | + |F 3 | ' U 

where the [ii are the magnifications and the Fi the ob- 
served fluxes, both with signs indicating the image pari- 
ties. The parities can be determined unambiguously be- 
cause in any triplet of adjacent images, the two outer 
images have the same parity while the middle image has 
the opposite parity (see Schneider et al. 1992; Petters et 
al. 2001). The ideal cusp relation has the form R cusp = 0. 

Note that we have defined i? CU sp to be non-negative. 
Several recent studies (Schechter & Wambsganss 2002; 
Keeton 2003; Kochanek & Dalai 2003) have pointed out 
that small-scale structure tends to suppress negative- 
parity images more often than it amplifies positive-parity 
images, while global perturbations generally do not dis- 
tinguish between images with different parities. In an 
ensemble of lenses with flux ratio anomalies, skewness in 
the signed i? C usp distribution may therefore distinguish 
local from global perturbations. However, the statistical 
nature of this argument precludes its use in individual 
lenses. Since we seek a method of identifying anoma- 
lies in individual lenses, we consider only the unsigned 
quantity. 

We first study the cusp relation analytically using two 
families of lens potentials where it is possible to obtain 
exact solutions of the lens equation. In one family, the 
galaxy is assumed to be spherical but is allowed to have a 
general power law surface density profile £ cx r a ~ 2 and to 
have an external shear 7. In the other family, the galaxy 
is assumed to have an "isothermal" profile £ a r _1 but 
is allowed to have a complex angular structure, including 
shear; we specifically consider an ellipsoidal galaxy per- 
turbed by multipole density fluctuations. Appendix [B] 
describes the two families of lens potentials in detail and 
gives solutions for the positions and magnifications of im- 
ages corresponding to sources on a symmetry axis of the 
lens potential. 

Figure [3] shows i? CU sp versus the opening angle 9 and 
separation d of an image triplet, for various potentials 
with different radial profiles, ellipticities, and shears. In 
general, i? C usp is small when 9 and d are small (indicat- 
ing that the source is very near a cusp), and grows as 
9 and d grow (indicating that the source is moving far- 
ther from the cusp). The analytic results allow us to 
understand how departures from the ideal cusp relation 
depend on properties of the lens potential. We see that 
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Fig. 3. — The cusp relation residual i? CUS p as a function of the 
opening angle 8 and the separation d of an image triplet, plotted 
for various lens potentials using the analytic solutions to the lens 
equation for sources on the major axis of the potential. In panel 
(c), 7 > (7 < 0) represents a shear aligned with (orthogonal to) 
the major axis of the galaxy. 



radical changes in the radial profile of the lens poten- 
tial — from a — 1 (isothermal) to a = (point mass) — 
have a negligible effect on the cusp relation. By contrast, 
moderate changes in the ellipticity and shear can affect 
the cusp relation by tens of percent. The fact that the 
cusp relation is quite sensitive to ellipticity, moderately 
sensitive to shear, and not very sensitive to the radial 
profile makes sense: reasonable changes in the angular 
structure of the potential (e and 7) can affect nearby 
images quite differently, while reasonable changes in the 
radial profile cannot. Incidentally, we note that when 
considering fixed ellipticity and shear amplitudes, i? C usp 
can be larger when the two are orthogonal than when 
they are aligned. 

The effects of multipole density perturbations are 
shown in Figure for lens potentials with an "isother- 
mal" (a = 1) radial profile. Multipole modes with m = 3 
or 4 and amplitudes of a few percent are common in the 
isophotes of observed early-type galaxies (Bender et al. 
1989; Saglia et al. 1993; Rest et al. 2001) and in the iso- 
density contours of simulated galaxy merger remnants 
(Heyl, Hernquist, & Spergel 1994; Naab & Burkert 2003; 
Burkert & Naab 2003); in particular, m = 4 modes with 
amplitudes 04 > can represent small disk-like compo- 
nents in bulge-dominated galaxies, which are not unusual 
(Kelson et al. 2000; Tran et al. 2003). Such modes might 
have a significant effect on the magnifications of lensed 
images (Evans & Witt 2002; Moller et al. 2003). We find 
that m = 4 modes do not significantly increase i? C usp 
for cusp triplets with 9 < 90° when the source is on the 
major axis of the lens potential (Figure^). However, 
they can create remarkably large values of i? CU sp even 
for small 9 when the source is on the minor axis (Fig- 
ure^). At fixed amplitude, higher order modes pro- 
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Fig. 4. — Effects of multipole perturbations on the cusp rela- 
tion, for isothermal ellipsoid lens potentials with e = 0.3 and 7 = 0. 
The heavy dashed curves all show a reference case with no mul- 
tipole modes, (a) The effects of m = 4 density fluctuations with 
amplitude 04 = (±0.01, ±0.02, ±0.03) indicated by increasing line 
thickness. Solid (dotted) curves correspond to 04 > (04 < 0). 
The sources lie on the major axis of the lens potential, (b) Similar 
to (a), but for sources on the minor axis of the potential (c) The 
effects of multipole perturbations of different orders, all with am- 
plitude a m = 0.02. Lines of increasing thickness and alternating 
type indicate m = (6,8,10,12,14,16). The cross on each curve 
marks the point with 8 = 640° /m. The sources lie on the major 
axis of the lens potential. 



duce progressively larger values of i? C usp at smaller an- 
gles (Figure^];). 8 The position of the peak in the i? C usp 
curve for different values of to can be approximated as 
# P eak ~ 640°/m. This result describes our fiducial case 
with e = 0.3 and a m = 0.02; varying e and 04 has a small 
(<5%) effect on the position of the peak, but a large ef- 
fect on the amplitude of the peak. Thus, we can say 
as a rule of thumb that image triplets with angle 9 are 
significantly affected only by modes with m > 640°/$. 
We conclude that it is important to consider multipole 
effects in the cusp relation analysis. But as it is not 
clear that real galaxies have percent-level perturbations 
in modes beyond to 4, it is equally important to hold 
the perturbations to reasonable levels. 

So far we have studied only sources lying on a symme- 
try axis of the lens potential. For the more general case 
we turn to Monte Carlo simulations. We pick random 
source positions and solve the lens equation (using the 
algorithm and software by Keeton 2001b) to generate a 

d, and i? C us P for 



catalog of mock lenses. We compute 




H 100 200 300 
opening angle 6 (deg) 



i 2 

separation d/R ei 



8 For the high-order multipole modes we do not show sources 
on the minor axis of the potential, because on the minor axis the 
caustics often have complicated butterfly catastr ophes that need 
not satisfy the cusp relation (see Appendix IB B.2l . 



, for various values of the ellipticity 
and shear and the angle A0 7 between them. The points show 
results for Monte Carlo simulations of random source positions. 
The heavy curves show the analytic results for on-axis sources, 
where solid (dotted) curves indicate the major (minor) axis of the 
lens potential. Analytic results are available only for A6 7 = and 
90°. 



each triplet in each 4-image lens, and then plot i? C usp 
versus 9 or d for all triplets. Figure El shows sample re- 
sults for isothermal ellipsoid galaxies with shear. The 
most important result is that over the region of interest 
for the cusp relation (9 < 180° and d/R c \ n < 1.7) there 
is a firm upper envelope on the values of i? CU sp- 9 In fact 
there are two envelopes: one each for major and minor 
axis cusps. Moreover, in lenses with reflection symme- 
try the envelope corresponds to sources on the symme try 
axis. To understand this result, in Appendix IB B.ll we 
prove that i? C usp is a local maximum on the symmetry 
axis of an isothermal sphere plus shear. Messy algebra 
hinders a rigorous analysis of other potentials, but intu- 
ition and the Monte Carlo simulations suggest that the 
result is generally true. In other words, the analytic re- 
sults for on-axis sources provide a simple and important 
upper bound on R CU sp- 

To summarize, the ideal cusp relation breaks down for 
sources a small but finite distance from the cusp, but in a 
way that can be understood and quantified. The realistic 
cusp relation is mainly sensitive to the angular structure 
of the lens potential, not the radial profile. The impor- 
tant quantities are the ellipticity, shear, and strength of 
multipole density fluctuations. For the subset of cusps 
that possess a symmetry axis, sources on that axis pro- 
vide a strict upper bound on i? CU sp over the interesting 
range of 9 and d that can often be derived analytically. 

9 The break at d/ij c ; n ~ 1.7 is simple to understand. This sep- 
aration corresponds to an image triplet comprising an equilateral 
triangle inscribed within the Einstein ring. When the separation 
reaches this value the images are so spread out that they can no 
longer be associated with a cusp. 
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5. THE CUSP RELATION IN REALISTIC LENS POPULATIONS 

If we knew the ellipticity, shear, and multipole per- 
turbations for individual observed lenses, we could use 
the previous analysis to compute how much R cusp can 
deviate from zero for smooth potentials and then con- 
clude that larger values represent flux ratio anomalies. 
Unfortunately, the three key quantities are not directly 
observable. The ellipticity and multipole perturbations 
of the mass need not be the same as those of the light 
(e.g., Keeton, Kochanek & Falco 1998), and in any case 
the shear cannot be directly observed. The three quanti- 
ties could be constrained with lens models, but we seek to 
avoid explicit modeling to the extent possible. Instead, 
our approach is to adopt observationally-motivated pri- 
ors on the ellipticity, multipole perturbations, and shear, 
and use Monte Carlo simulations to obtain a sample of 
realistic lens potentials and derive probability distribu- 
tions for i?cusp- In this section we describe the priors 
( 4b 5.1fl and methods ( 4b 5.2|) for the simulations. 

5.1. Input distributions 

We consider only early-type galaxies, because they are 
expected to dominate the lensing optical depth due to 
their large average mass (e.g., Turner, Ostriker, & Gott 
1984; Fukugita & Turner 1991). Indeed, -80-90% of ob- 
served lens galaxies have properties consistent with being 
massive ellipticals (Keeton et al. 1998; Kochanek et al. 
2000; Rusin et al. 2003a). The distinction between ellip- 
ticals and spirals is important, because disk-dominated 
galaxies that are viewed close to edge-on can produce 
cusp configurations that deviate significantly from the 
cusp relation (Keeton & Kochanek 1998). Several of the 
lenses for which we identify flux ratio anomalies are con- 
firmed ellipticals, and none of them have properties sug- 
gesting that they are spirals (Impey et al. 1996; Burud et 
al. 1998; Fassnacht et al. 1999; Jackson et al. 2000; Inada 
et al. 2003; Sluse et al. 2003; Rusin et al. 2003a). 10 

We allow the simulated galaxies to have ellipticity and 
also octopole (m = 4) perturbations, with distributions 
drawn from observations of isophote shapes in early-type 
galaxies. Even if the shapes of the mass and light dis- 
tributions are not correlated on a case-by-case basis, it 
seems likely that their distributions are similar (see Rusin 
& Tegmark 2001 for a discussion). Indeed, the distribu- 
tion of isodensity contour shapes in simulated merger 
remnants is very similar to the observed distribution 
of isophote shapes (Heyl et al. 1994; Naab & Burkert 
2003; Burkert & Naab 2003). Multipole perturbations 
beyond to > 5 have generally not been reported, but it 
is likely that they must have relatively low amplitudes 
to be compatible with observations. Lower-order m = 3 
modes have been reported with amplitudes comparable 
to to = 4 modes (e.g., Rest et al. 2001), but we do not 
consider them here because they are not reported in the 
samples we use, and because at fixed amplitude higher- 
order modes produce larger deviations in the cusp rela- 
tion (see Figure 0J). Our approach is formally equivalent 
to studies that explicitly include a disk-like mass com- 

10 Despite a suggestion by Moller et al. (2003) that the lens 
galaxy in B2045+265 might be a spiral, its structural and dynam- 
ical properties are fully consistent with being an elliptical (Rusin 
et al. 2003a), and no disk-like structure is evident in Hubble Space 
Telescope images (C. Kochanek, private communication). 
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Fig. 6. — The main panel shows the e and values for the 
galaxies in the Bender et al. (1989) and Saglia et al. (1993) sam- 
ples. Galaxies with 04 > (124 < 0) have disky (boxy) isophotes. 
The right panel shows the histograms of 04. The top panel shows 
histograms of e for these samples and also the j0rgensen et al. 
(1995) sample. 



ponent (e.g., Moller et al. 2003), since small disks can 
be treated as to = 4 multipole perturbations. We use 
only isothermal galaxies, because the radial profile of the 
lens galaxy does not significantly affect the cusp relation. 
The galaxy mass is unimportant, because it simply sets 
a length scale (the Einstein radius) that can be scaled 
out by using the dimcnsionless separation d/R c - m . 

We use ellipticity and octopole distributions from three 
different samples. Jorgensen et al. (1995) report elliptic- 
ities for 379 E and SO galaxies in 11 clusters, including 
Coma. Their ellipticity distribution has mean e = 0.31 
and dispersion o~ e = 0.18. Since the Jorgensen et al. sam- 
ple does not include octopole amplitudes, we consider 
two smaller samples that do. Bender et al. (1989) report 
cllipticities and octopoles for 87 nearby, bright ellipti- 
cal galaxies. Their ellipticity distribution has e = 0.28 
and o~ e — 0.15, while their octopole distribution has 
mean a 4 = 0.003 and dispersion a ai — 0.011. Finally, 
Saglia et al. (1993) report ellipticities and octopole am- 
plitudes for 54 ellipticals in Coma, with e = 0.30 and 
<j e = 0.16, and 04 = 0.014 and a a4 — 0.015. Compared 
to the Bender et al. sample, the Saglia et al. sample has 
a higher incidence of galaxies with strong disky pertur- 
bations (04 > 0). Considering all three samples allows 
us to examine whether our conclusions depend system- 
atically on the input (e, 04) distributions (although we 
note that since the Jorgensen et al. and Saglia et al. 
samples both include galaxies in Coma they are not fully 
independent). Figure shows the different samples and 
suggests that e and a 4 are correlated such that highly 
elliptical galaxies tend to have significant disky pertur- 
bations. Our analysis includes this correlation explicitly 
by using the observed joint distribution of e and a 4 . 
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For the shear amplitude we adopt a lognormal distri- 
bution with median 7 = 0.05 and dispersion cr 7 = 0.2 
dex. This is consistent with the distribution of shears 
expected from the environments of early- type galaxies, 
as estimated from .ZV-body and semi-analytic simula- 
tions of galaxy formation (Holder & Schechter 2003). 
It is broadly consistent with the empirical distribution 
of shears required to fit observed lenses, when selection 
biases related to the lensing cross section and magnifica- 
tion bias are taken into account (see Holder & Schechter 
2003). The mean shear is also consistent with the typ- 
ical value needed to explain misalignments between the 
light and mass in observed lenses (Kochanek 2002). We 
assume random shear orientations. 

5.2. Simulation methods 

We use each input distribution to run a large Monte 
Carlo simulation containing ~ 10 6 4-image lenses. With 
the J0rgensen et al. sample we draw 2000 ellipticities 
from the observed ellipticity distribution and give each a 
random shear. With the Bender et al. and Saglia et al. 
samples we use only the observed (e,a4) pairs to make 
sure we include the apparent correlation between the two 
quantities, but we use each pair with 100 different ran- 
dom shears; thus we consider 8700 and 5400 lens poten- 
tials for the Bender et al. and Saglia et al. samples, re- 
spectively, but we need to remember that these represent 
only 87 and 54 different ellipticity and octopole measure- 
ments. For each potential, we pick random sources with 
a uniform density of ~ 10 3 R~ in in the source plane and 
solve the lens equation using the algorithm and software 
by Keeton (2001b). We have verified that our results 
are not sensitive to the number of shears and density of 
sources used. 

To understand how our mock lenses compare to ob- 
served samples it is important to consider two selection 
effects. First, the cross section for 4-image lenses is very 
sensitive to ellipticity and shear, but our uniform sam- 
pling of the source plane ensures that each lens potential 
is automatically weighted by the correct cross section. 
Second, magnification bias can favor lenses with higher 
amplifications. While this effect is important when com- 
paring 4-image lenses to 2-image lenses (e.g., Keeton et 
al. 1997; Rusin & Tegmark 2001), it is much less im- 
portant when comparing different 4-image lenses against 
each other. If anything, it would favor the sources very 
near a cusp or fold that yield extremely magnified lenses 
that best satisfy the cusp/fold relations, giving more 
weight to lenses with smaller deviations from the ideal 
relations. We therefore neglect magnification bias and 
believe that this is a conservative approach. 

We compute 9, d, and i? C usp for each image triplet in 
each mock 4-image lens, and use this ensemble to deter- 
mine the conditional probability distributions p(R cusp \9) 1 
p(Rcusp\d), and p(i? cusp |<i, 0) — the probability of having 
a particular value of i? C usp given 9, d, or both. For exam- 
ple, Figure [7| shows curves of constant conditional proba- 
bility p(i? C usp| 9) and p(R cusp \d) versus 9 and d. This fig- 
ure is basically a modified version of Figure El where we 
have averaged over appropriate ellipticity, octopole, and 
shear distributions. It is interpreted as saying that 68% 
of triplets in our sample of mock lenses lie in the region 
between the solid curves, 99% lie between the dashed 
curves, and so forth. To the extent that our simula- 



tions encompass the range of ellipticities, octopoles, and 
shears in real populations of early-type galaxies, we can 
conclude that any points lying outside the contours rep- 
resent flux ratios that are inconsistent with smooth lens 
potentials. 

6. APPLICATION TO OBSERVED LENSES 

We can now use our theoretical analysis to evaluate ob- 
served lenses, seeking to identify systems that violate the 
cusp relation and therefore have anomalous flux ratios. 
We first summarize the data f M6 6.1|) and then present 
our results f H6 6.2fl . 

6.1. Data 

The data for the nineteen published 4-image lenses are 
given in Table IB II Only five of the lenses are thought 
to have cusp configurations, but we can still apply the 
cusp relation analysis to all of them to see what we learn. 
In the table we list the four different image triplets for 
each lens, with the opening angle 9, the image separa- 
tion d, and the observed value of i? C usp for each. If the 
lens galaxy position is known, the angle 9 is fully de- 
termined by the data; if not, we estimate 9 using the 
galaxy position estimated from lens models, which is a 
fairly model-independent prediction. The separation d 
is determined directly from the data. To normalize it 
we need the Einstein radius i? e in, which must be de- 
rived from a lens model but is insensitive to the assumed 
model; different models generally yield the same Einstein 
radius with systematic uncertainties of just a few percent 
(Cohn et al. 2001; Rusin et al. 2003b). 

When we measure i? C usp we need to consider system- 
atic uncertainties due to effects like source variability and 
the lens time delay, scatter broadening (at radio wave- 
lengths), and differential extinction by patchy dust in 
the lens galaxy (at optical wavelengths). (Extinction by 
dust in our own Galaxy does not affect the flux ratios, 
because it affects all images equally.) Dalai & Kochanek 
(2002) advocate adopting a fiducial estimate of 10% un- 
certainties in the flux ratios to account for these effects, 
but this is likely to be quite conservative. For most lenses 
the uncertainties are irrelevant because the measured val- 
ues of -Rc usp h e well within the expected distribution, so 
for Table IbD we use 10% flux uncertainties for simplic- 
ity. We want to be more careful about the error budgets 
for lenses suspected of having flux ratio anomalies, so we 
discuss them individually in t he next section. 

For comparison, Table IB II also gives values for R cusp 
predicted by standard lens models consisting of an 
isothermal ellipsoid with an external shear. 11 Only 
the image positions (not the flux ratios) were used as 
constraints. In MG 0414+0534, RX J0911+0551, and 
B1608+656 the lens models also include the perturba- 
tive effects of an observed satellite galaxy near the main 
lens galaxy. We stress that the lens models are not actu- 
ally used in seeking flux ratio anomalies (other than for 
estimating i? c in, as discussed above). They are included 
only as a general indication of what to expect for i? C usp 
from smooth lens models for these systems. 

11 The only exception is B1555+375, where the ellipsoid plus 
shear model is somewhat ambiguous and we use models with a 
slightly different parameterization of the quadrupole moment of 
the lens potential; see Turner et al. (C. Turner, C. R. Keeton, & 
C. S. Kochanek, in prep.) for technical details. 
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Fig. 7. — The curves show contours of constant conditional probability p(R C ubp\8) (top) and p(i? C uspM) (bottom), from Monte Carlo 
simulations using the J0rgcnsen et al. input data (left) or the Bender et al. input data (right). The contours are drawn at the 68% 
(solid), 95% (long-dashed), 99% (dashed), and 99.9% (dotted) confidence levels. The points show the observed values of i? CU sp for the 
known 4-image lenses. B0712+472 appears twice; the lower value of -Rcusp corresponds to the radio data, while the higher value (labeled) 
corresponds to the optical/near-IR data (see text). 



6.2. Results 

Figure shows the observed values of i? CU sp super- 
posed on the predicted confidence contours. The five 
cusp lenses can be identified as the points with 9 < 80° 
(B0712+472 appears twice: once for radio data, and once 
for optical/near-IR data). Most of the observed lenses 
lie within the predicted confidence region, so according 
to this analysis they are not obviously inconsistent with 
smooth lens potentials. We note that the predicted con- 
fidence contours are very similar for the J0rgensen et al. 
and Bender et al. input data (and also for the Saglia et al. 
input data, not shown). The main difference is that the 
presence of octopole perturbations in the Bender et al. 
input data causes the confidence contours to stretch to 
higher values of i? C usp at d/R e { n < 1.7, which is what we 
expect from the theoretical analysis in 21 (The contours 
for the Bender et al. input data are somewhat noisy due 
to the relatively small number of ellipticity and octopole 
measurements.) Thus, contrary to the claim by Moller 
et al. (2003), we find that adding (properly- normalized) 
disky components to elliptical galaxies does not have an 
enormous effect on the cusp relation. We shall explain 
below which of our conclusions are or are not affected 
by the presence of octopole terms, or more generally by 
changes in the input data. 

Several of the lenses are obvious outliers. The cusp 
lens B2045+265 lies outside all contours. The cusp lens 
B0712+472 lies outside all p(R C us P \9) contours and ei- 
ther outside or just inside the 99.9% confidence contour 
for p(R cusp \d), depending on the input data. The cusp 
lenses 1RXS J1131-1231 and RX J0911+0551 stand out 
relative to p(R C usp\d) for the J0rgensen et al. input data 
but not for the Bender et al. (or Saglia et al.) input data. 
(Incidentally, RX J0911+0551 and B2045+265 are also 



responsible for the p(R CU sp\9) outliers at 9 — 290° and 
9 = 325°, respectively.) The fifth cusp lens B1422+231 
does not stand out in this analysis. Finally, the lens 
SDSS J0924+0219 is an outlier with respect top(i? cusp |6>) 
even though it is not a cusp configuration. 

The joint conditional probability distribution 
p(R cuap \d, 9) provides an even more powerful way 
to identify outliers. Figure |H1 compares the cumulative 
probability -P m od(> Rcusp\d, 9) that smooth potentials 
produce R cusp larger than some value versus the cumu- 
lative probability P bs(< -Rcusp) that the measurement 
of i? C usp is smaller than some value, for the six lenses 
just mentioned. The measured value of i? C usp is com- 
patible with smooth potentials only if the curves have 
a significant overlap. B2045+265, 1RXS J1131-1231, 
and SDSS J0924+0219, and B0712+472(optical) are 
clear outliers; B0712+472(radio) and RX J0911+0551 
are marginal outliers; and B1422+231 is the only case 
where the observed and predicted distributions are 
clearly compatible. We now discuss each of these lenses 
individually. 

6.2.1. B2045+265 

Fassnacht et al. (1999) give eight measurements of the 
radio fluxes for B2045+265 at 1.4, 5, 8.5, and 15 GHz, 
from different radio arrays with different resolutions. 
The mean value and scatter in R cusp is 0.516 ± 0.018; 
the scatter is only slightly larger than the uncertainty 
that would be inferred from the quoted flux errors. The 
fact that the i? CU sp values from diverse radio datasets 
are consistent within the errors argues against any sig- 
nificant non-gravitational effects (e.g., scattering). Also, 
for a cusp triplet the time delays are expected to be very 
short — predicted to be <6 hours for B2045+265, and 
similarly short for the other cusp lenses — so they should 
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discussed individually. The rising curves show the probability 
^obs(<^cusp) that the measurement of i? C usp is smaller than some 
value; for B0712+472 the solid (dashed) curves denote the radio 
(optical/near-IR) measurements. The falling curves show the prob- 
ability P mo d(> Rcusp\d,9) that smooth lens potentials produce 
^icusp larger than some value, with different line types denoting 
different input data for the Monte Carlo simulations. 



have no effect on the measured value of i? C us P - We there- 
fore believe that ±0.018 represents a reasonable estimate 
of the uncertainty. Koopmans et al. (2003a) present 41 
measurements of B2045+265 at 5 GHz. Although they 
observe variability that they attribute to scintillation, 
they find a mean and scatter in R cusp of 0.501 ± 0.035 in 
excellent agreement with the value from the Fassnacht et 
al. (1999) data. 

The Cf A/ Arizona Space Telescope Lens Survey (CAS- 
TLES; C. Kochanek et al., private communication) 12 
provides data at optical and near-IR wavelengths from 
Hubble Space Telescope imaging. Their data for 



B2045+265 yield R, 



cusp 



0.501 ± 0.037 in V-band, 
0.531 ± 0.035 in I-band, and 0.502 ± 0.015 in H-band. 
The colors of the images, and the fact that R cusp re- 
mains constant over a factor of three in optical/near-IR 
wavelength, indicate that there is little or no differential 
extinction between the images. The weighted average 
of the optical/near-IR data yields -R CU sp = 0.506 ±0.013; 
the excellent agreement with the radio data suggests that 
the measured value of -R CU sp is robust and independent 
of wavelength, and that the small inferred uncertainties 
on -Rcusp are realistic. The weighted average of all mea- 
surements is i?cusp = 0.509 ± 0.010. 

Figures and |Hl show that the existence of a flux ra- 
tio anomaly in B2045±265 is beyond doubt. Image B is 
simply much too faint to be consistent with smooth lens 
potentials, no matter which input ellipticity and octopole 
distributions are used. Attempting to explain the value 
of -Rcusp with multipolc perturbations would require a 
significant amplitude in a mode with m 
urelU). 



6.2.2. 



16 (see Fig- 



B071 2+472 

Jackson et al. (1998) give three different measurements 
of the radio fluxes for B0712+472 at 5 GHz and one 

12 See http://cfa-www.harvard.edu/castles 



measurement at 15 GHz. The mean and scatter in the 
value of -Rcusp from the four datasets is 0.261 ± 0.031. 
The data from 41 measurements at 5 GHz by Koopmans 



et al. (2003a) yield R c 



0.255 ± 0.030, in excellent 



agreement with the Jackson et al. (1998) value. The 
weighted average of these measurements is 0.258 ±0.022. 

The optical and near-IR data from CASTLES (also see 
Jackson, Xanthopoulos, & Browne 2000) yield -R CU sp = 
0.619 ± 0.050 in V-band, 0.572 ± 0.147 in I-band, and 
0.473 ± 0.092 in H-band. The decline in R cusp with 
wavelength suggests that there might be some differen- 
tial extinction between the images, but the evidence is 
weak because the three measurements are formally con- 
sistent within the errors. The weighted average of the 
optical/near-IR measurements is 0.585 ± 0.042. 

The difference between the radio and optical results is 
very interesting. At radio wavelengths the overlap be- 
tween the observed and predicted probability curves in 
Figure |H1 is small but non-negligible: the curves overlap 
at 2.1-3.7% depending on the input data used in the 
Monte Carlo simulations. Thus, there is evidence for a 
radio flux ratio anomaly, but only at the 96-98% confi- 
dence level. At optical/near-IR wavelengths, by contrast, 
there is no overlap between the observed and predicted 
probability curves, and hence evidence for an optical flux 
ratio anomaly at high confidence. Both conclusions are 
unaffected by octopole perturbations in the lens poten- 
tial, and more generally by changes to the input data in 
the Monte Carlo simulations. The difference between the 
radio and optical results could indicate that the optical 
flux ratio anomaly is caused by a star (microlensing) or 
some other object with a characteristic size smaller than 
a typical dark matter subhalo. That possibility makes 
B0712+472 a promising system for optical monitoring to 
look for variability that would indicate microlensing. 

6.2.3. 1RXS J1131-1231 

Sluse et al. (2003) present three measurements of the 
optical flux ratios of 1RXS Jl 131 — 1231. Observations 



from 2 May 2002 yield R c 



0.350 ±0.021 in V-band, 



while observations from 18 December 2002 yield -R CU sp = 
0.353 ± 0.031 in V-band and 0.367 ± 0.031 in R-band. 
It is interesting that the total flux of the system varied 
by 0.29 ± 0.04 mag between May and December, yet the 
flux ratios and -R CU sp values are essentially identical. The 
likely explanation is that the source varied over the 7- 
month time scale, but the short time delays between the 
bright images A, B, C (predicted to be <1 day) kept 
the flux ratios essential constant. The fact that -R CU sp is 
the same at different epochs and in different passbands 
indicates that there are no large systematic uncertainties 
due to the time delays or electromagnetic effects. 

The weighted average of the -R CU sp values is 0.355 ± 
0.015. Figure |H1 shows that this value lies well outside 
the predicted distributions, implying a strong flux ra- 
tio anomaly, and that this conclusion is insensitive to 
changes in the input data in the Monte Carlo simula- 
tions. 



6.2.4. RX J0911+0551 

RX J0911±0551 is the only cusp lens that shows sig- 
nificant evidence for differential extinction between the 
images. In the CASTLES data (also see Burud et al. 
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Fig. 9.— Results for RX J0911+0551 for the "normal shear 
and "cluster shear" cases, using the Bender et al. input data. 



1998), image Ai has colors V - H = 1.24 ± 0.04 and 
7--ff = 0.79±0.04 while images A 2 and A 3 both have 
colors V - H = 1.54 ± 0.05 and I — H = 1.01 ± 0.04. At- 
tributing the color difference to dust in the lens galaxy, 
we estimate a differential extinction between Ai and the 
A 2 /A 3 pair of (0.44,0.32,0.12) mag in (V,I,H) (fol- 
lowing the analysis of Falco et al. 1999, using a red- 
shifted Ry = 3.1 extinction curve from Cardelli, Clay- 
ton & Mathis 1989). 13 Correcting for the reddening 
changes R cusp from 0.226 for the raw H-band flux ratios 
to i?cusp = 0.192 ± 0.011 for the de-reddened data. 

Figure [S] compares this value with the predictions of 
smooth lens potentials. Perhaps the most interesting 
point is that RX J0911+0551 is the only system where 
octopole perturbations affect the judgment about a flux 
ratio anomaly. The Monte Carlo simulations that include 
octopole modes (based on the Bender et al. or Saglia et 
al. input data) have a long tail to high -R CU sp values that 
is absent from simulations lacking such modes (based on 
the J0rgensen et al. input data, or re-running the Ben- 
der et al. or Saglia et al. data with the octopole terms 
removed [not shown]). The octopole modes apparently 
allow smooth lens potentials to create a small but sig- 
nificant number of image triplets with the same angle 
9 and separation d as RX J0911+0551 that have rel- 
atively high values of -R CU sp- The reason we see such 
an effect here, but not in the other lenses we study, 
is that RX J0911+0551 is the only system where the 
source lies near a cusp on the minor axis of the galaxy. 
As we showed in 21 octopole perturbations mainly af- 
fect sources near a minor-axis cusp. The octopole modes 
cause the observed and predicted probability curves in 
Figure |H1 overlap at 1.7% for the Bender et al. input data 
or 8.6% for the Saglia et al. input data. We caution 
against overinterpreting these numbers, however. Be- 
cause the Bender et al. and Saglia et al. samples include 
just 87 and 54 ellipticity and octopole measurements, 
respectively, we may worry about small-number statis- 
tics and sample variance. Indeed, the difference between 
the i?cusp distributions predicted by these two samples 
suggests that a robust analysis of the cusp relation for 
minor-axis lenses like RX J0911+0551 will require larger 

13 Moller et al. (2003) claim that the extinction correction for 
RX J0911+0551 is highly uncertain, but we find that not to be the 
case. 



samples with better determinations of the ellipticity and 
multipole distributions. 

The analysis of RX J091 1+0551 is actually even more 
complicated, because the lens galaxy lies in a cluster en- 
vironment that contributes a large shear 7 ~ 0.3 to the 
lens potential (Kneib, Cohen, & Hjorth 2000). Since this 
shear is almost 4tr above the median of our assumed dis- 
tribution, it is not likely to be well represented in our 
fiducial Monte Carlo simulations. To examine the effects 
of such a large shear, we use a new Monte Carlo simu- 
lation with the shear amplitude fixed at 7 = 0.3. The 
orientation of the shear is still random, and we adopt 
the Bender et al. (1989) input data for the ellipticity 
and octopole. Figure [§] shows that lens potentials with 
a large shear produce a narrower range of R CUS p values 
than lens potentials with the fiducial shear distribution. 
The likely explanation is that an octopole term has the 
most effect when aligned with the quadrupole moment of 
the potential, which cannot happen when the quadrupole 
is dominated by a large, randomly oriented shear. As a 
result, the -R CU sp distribution predicted when the shear is 
large no longer overlaps with the observed value of -R CU sp- 

We believe that the "cluster shear" case pro- 
vides strong evidence for a flux ratio anomaly in 
RX J0911+0551. However, it will be important to exam- 
ine more carefully whether multipole modes can compro- 
mise this conclusion. The most straightforward approach 
will be direct lens modeling. It will definitely be possi- 
ble to fit the current data exactly by allowing enough 
multipole modes in the models (see Evans & Witt 2002 
for examples). The crucial test would be to obtain deep 
near-IR imaging to try to find an Einstein ring image of 
the quasar host galaxy (e.g., Impey et al. 1998; Keeton 
et al. 2000; Kochanek et al. 2001). The additional con- 
straints from the Einstein ring would greatly restrict the 
space of allowed models and determine whether the mul- 
tipole modes required by the flux ratios are acceptable 
or not (see Kochanek & Dalai 2003). 

If the putative flux ratio anomaly in RX J0911+0551 
is confirmed, its interpretation may be ambiguous. With 
only broad-band optical and near-IR flux ratios it is im- 
possible to determine whether the small-scale structure 
implied by the anomaly corresponds to micro-lensing or 
milli-lcnsing — or to one or more dwarf galaxies, which 
tend to be abundant in cluster environments (Trentham 
1997, Smith, Driver, & Phillipps 1997; Driver, Couch, & 
Phillipps 1998; Zabludoff & Mulchaey 2000). In fact, the 
CASTLES images do reveal a faint satellite galaxy lying 
l'.'O away from the main lens galaxy (although not close 
to any of the lensed images). Analysis of lens models 
that explicitly include the satellite suggest that it can- 
not explain the observed flux ratios; as seen in Table |B~T1 
such models still predict R CU s P < 0.01 for the cusp triplet. 
Still, we cannot presently rule out the possibility that the 
anomaly is caused by an as- yet undetected but otherwise 
unremarkable dwarf galaxy. Perhaps the best prospect 
for determining whether the anomaly is caused by a star, 
a dark matter clump, or a dwarf galaxy would be high- 
resolution resolved spectroscopy of the three images (see 
Moustakas & Metcalf 2003). 

6.2.5. B1422+231 

Patnaik & Narasimha (2001) report 61 measurements 
of the radio fluxes for B1422+231 at 8.4 and 15 GHz. The 
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mean and scatter for i? C usp from all of the measurements 
is 0.179 ± 0.006. The error on i? C us P obtained directly 
from the individual flux errors is comparable to or smaller 
than the scatter. Koopmans et al. (2003a) report 41 
measurements at 5 GHz that yield B, cusp = 0.187 ±0.004, 
in reasonable agreement with the value from the Patnaik 
& Narasimha (2001) data. 

The optical and near-IR data from CASTLES yield 
R cusp = 0.223 ±0.055 in V-band, 0.222 ±0.038 in I-band, 
and 0.175 ±0.015 in H-band. The change from V/I to H 
is comparable to the uncertainties. The weighted average 
of the optical/near-IR results is i? C us P = 0.184 ± 0.014, 
in remarkably good agreement with the radio results. 

Figures [7| and |H1 show that B1422+231 is not an outlier 
in the cusp relation according to our analysis. This re- 
sult is interesting in light of detailed modeling that shows 
B1422+231 to have a flux ratio anomaly (Mao & Schnei- 
der 1998; Bradac et al. 2002; Metcalf & Zhao 2002). The 
reason for the difference is that we are focusing on the 
maximum allowed values of i? C usp for given values of 9 
and d. We do not consider whether the potentials that 
can produce those values are at all consistent with the 
rest of the lens data. In other words, we are throwing 
away reliable information in our attempt to be as general 
and robust as possible. Metcalf & Zhao (2002) demon- 
strate that imposing the position constraints via detailed 
lens modeling provides a more sophisticated analysis that 
has more power to identify flux ratio anomalies. But be- 
cause it relies on specific modeling, we believe that such 
model fitting is less robust and generic than the cusp 
relation analysis. 

6.2.6. SDSS J0924 +0219 

The lens SDSS J0924+0219 is not a typical cusp config- 
uration but rather more like a cross configuration. Never- 
theless, it does have what seems to be an enormous flux 
ratio anomaly. The anomaly appears as such a strong 
suppression of one of the negative-parity images that the 
system was initially thought to have just three images, 
and the fourth image (D) was identified only after sub- 
traction of the three main images and the lens galaxy 
(Inada et al. 2003). Because image D is a factor of 14 
fainter than image A, the triplet ABD (with 9 — 156° 
and d/R em = 2.1) has i? C us P values of 0.916 ± 0.009 in 
5-band, 0.903 ± 0.005 in r-band, and 0.894 ± 0.005 in 
i-band. The good agreement between the different pass- 
bands suggests that there is little differential extinction 
between the images. The weighted average of the mea- 
surements is i?cus P = 0.901 ± 0.003. 

The fact that SDSS J0924+0219 is not a tight cusp 
configuration means that smooth potentials can produce 
a fairly broad distribution of i? C usp values (see Figure |HJ. 
Nevertheless, the observed value of i? C usp is so large that 
it still lies well outside the predicted distribution. In 
other words, the anomaly in SDSS J0924±0219 is so 
strong that it is identified by the cusp relation analy- 
sis even though the lens is not a proper cusp configura- 
tion. This conclusion is insensitive to changes in the in- 
put data of the Monte Carlo simulations, and is therefore 
robust provided that the lens galaxy is early-type. Be- 
cause the smallest image separation is still relatively large 
(d/i? c in = 2, 9 = 156°), a spiral galaxy might in prin- 
ciple have structure on scales appropriate to explain the 
anomaly. However, the strength of the anomaly would 



require significant structure, which would probably affect 
not just image D but all four images; and there is no cur- 
rent evidence that the lens is a spiral (Inada et al. 2003). 
In any case, the flux ratio anomaly in SDSS J0924+0219 
indicates that there is some strong and interesting struc- 
ture in the lens potential that is generally incompatible 
with the known properties of the luminous components 
of early- type galaxies. 

6.3. Comments 

To conclude this section, we remark that it is surprising 
to see how large i? C usp can be for realistic lenses, even in 
the absence of small-scale structure. Image triplets must 
be quite tight (small 9 and d) in order to be guaranteed 
of satisfying the cusp relation reasonably well. For ex- 
ample, to satisfy the relation with i? cusp < 0.1 at 99% 
confidence, the opening angle must be 9 < 30°. The 
parameter space that gives rise to such image configu- 
rations is quite small; in our Monte Carlo simulations, 
only ~0.1% of 4-image lenses have cusp triplets that are 
this tight. While this estimate omits magnification bias 
(which will favor tight configurations), it does suggest 
that cusp image configurations that are close to ideal are 
likely to be rare. The cusp relation can still be a valuable 
tool for identifying flux ratio anomalies, but it must be 
used with some care. 

We have identified flux ratio anomalies that cannot 
satisfy the cusp relation for any (reasonable) combina- 
tion of shear, ellipticity, and octopole modes. While we 
do not believe that significant high-order multipole per- 
turbations are common in real early- type galaxies, it is 
nevertheless interesting to consider what modes would 
be required to explain the anomalies. Using the ana- 
lytic results from ^1 (Figure 0] in particular) , we esti- 
mate that B2045±265 would require 1.5% density fluctu- 
ations in a mode with m > 16. The B0712+472 optical 
anomaly would require 4.5% fluctuations with m > 8. 
1RXS J1131— 1231 would require 1.5% fluctuations with 
m > 10. SDSS J0924±0219, being a cross rather than a 
cusp configuration, does not require a particularly high- 
order mode (m > 4), but it would require a perturbation 
amplitude of at least 12%. The estimated amplitudes 
are based on an assumed ellipticity e = 0.3, and they 
would increase or decrease as the ellipticity is increased 
or decreased. Note that these estimates are consistent 
with our rule of thumb from 2] that a large i? C usp value 
in a triplet with angle 9 requires multipole modes with 
m > 640°/6'. Thus, the most general statement about 
the lens potential that can be derived from a flux ratio 
anomaly is that there must be significant structure on 
the scale of the separation between the images. 

7. CONCLUSIONS 

The images associated with an ideal cusp catastrophe 
satisfy universal position and magnification relations, the 
most interesting of which says that the signed magnifi- 
cations of the three images should sum to zero. (The 
other relations have less practical value because they in- 
volve unobservable quantities.) A violation of this rela- 
tion indicates that the catastrophe is not an ideal cusp 
— that the lens potential contains terms beyond fourth 
order (third order in the lens equation; see eq. A 
significant violation of the cusp relation can imply that 
an observed 4-image lens has significant structure in the 
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lens potential on scales comparable to or smaller than 
the separation between the images. 

An important caveat is that the ideal cusp relation 
is not expected to hold exactly for real lenses, because 
real lens galaxies are not ideal cusps. It is therefore cru- 
cial to understand the degree to which simple aspects 
of real lens potentials (the radial density profile, ellip- 
ticity, multipole density fluctuations, and tidal shear) 
can cause deviations from the ideal cusp relation. Af- 
ter showing that the radial profile has little effect on 
the analysis, we have adopted observationally- motivated 
distributions for the shear, ellipticity, and octopole per- 
turbations and derived the resulting distribution for the 
quantity i? cusp = (|F\ + F 2 + F 3 |)/(|iq| + \F 2 \ + \F 3 \), 
where the Fi are the signed fluxes for a triplet of images. 
The distribution describes the range of R cusp values that 
can be produced by lens potentials containing realistic 
shears, ellipticities, and octopole modes. If an observed 
4-image lens has an i? CU sp value lying outside this range, 
we can conclude that its lens potential must have some 
kind of structure on scales comparable to or smaller than 
the image separation that is not represented in the lu- 
minous properties of early-type galaxies or the external 
shears expected from their typical environments. This 
structure might be large, low-order multipole moments, 
significant power in moments higher than octopole or- 
der, or some kind of compact structure. Analysis of the 
cusp relation therefore provides a new way to search for 
cases of milli-lensing, micro-lensing, or other interesting 
small-scale structure, which is an attractive alternative 
to global lens modeling because it relies on a local anal- 
ysis of lensing near a cusp to probe localized structure in 
the lens potential. 

When we examine the nineteen known 4-image lenses, 
we find evidence for flux ratio anomalies (violations of the 
cusp relation) in four of the five lenses with cusp image 
configurations, plus one other lens. First, B2045+265 
has a very strong flux ratio anomaly at both radio 
and optical/near-IR wavelengths, which if attributed to 
multipole density fluctuations would require significant 
structure in modes with m > 16. Second, B0712+472 has 
an anomaly that is very strong at optical/near-IR wave- 
lengths but marginal at radio wavelengths, which would 
require structure in multipole modes with m > 8. The 
fact that the anomaly is much stronger at optical/near- 
IR wavelengths than at radio wavelengths might suggest 
that it is due to micro-lensing rather than milli-lensing, 
but further study is required to test that hypothesis. 
Third, 1RXS J1131— 1231 has a strong anomaly at op- 
tical wavelengths, which would require structure in mul- 
tipole modes with m > 10. Fourth, RX J0911+0551 
has a value of R CUS p that suggests a flux ratio anomaly 
but whose interpretation is somewhat complicated. As 
the only known system where the source lies near a cusp 
on the minor (rather than major) axis of the lens po- 
tential, RX J0911+0551 is the only one where octopole 
modes can allow simple lens potentials to produce rel- 
atively high values of i? C us P . Strong conclusions about 
the putative anomaly will require better knowledge of 
the distribution of octopole amplitudes in real galaxies. 
Still, we believe that the evidence for a flux ratio anomaly 
in this system is good, especially when the fact that the 
lens lies in a cluster environment is taken into account. 
Finally, SDSS J0924+0219 is an intriguing system with 



a flux ratio anomaly so strong that it is identified by the 
cusp relation analysis even though the lens does not have 
a cusp image configuration. 

Interestingly, we find that the cusp lens B1422+231 
does not obviously violate the cusp relation, even though 
it is known to have a flux ratio anomaly from detailed 
lens modeling (Mao & Schneider 1998; Bradac et al. 2002; 
Metcalf & Zhao 2002). This illustrates an important 
point: in analyzing the cusp relation we have ignored 
constraints from the observed image positions, in order 
to be general and to avoid explicit modeling as much 
as possible. The idea is that an analysis based purely 
on magnification relations is the most robust and con- 
servative way to identify flux ratio anomalies. Adding 
constraints from the image positions yields an analysis 
that is more sophisticated and has more power to iden- 
tify flux ratio anomalies (see Metcalf & Zhao 2002 for 
a good example), but it requires detailed modeling and 
is therefore less generic. Furthermore, even when direct 
modeling suggests that a set of flux ratios cannot be fit 
by simple lens potentials, it may not reveal why that is 
the case. The cusp relation immediately pinpoints the 
cause of the failure. For these reasons, we believe that 
the cusp relation analysis is the best place to start when 
seeking to identify lenses with flux ratio anomalies (at 
least for lenses with cusp configurations) . If an observed 
lens violates the cusp relation then the anomaly is un- 
ambiguous and easy to understand. If it does not violate 
the cusp relation, it may still have an anomaly but may 
only be revealed by more sophisticated and less model- 
independent analyses. 

Although we have argued that flux ratio anomalies ex- 
ist and indicate a strong need for small-scale structure in 
lens galaxies, we cannot draw strong conclusions about 
what that structure must be. Analyses based on image 
positions and broad-band flux ratios can only give an up- 
per limit on the characteristic angular scale of the struc- 
ture implies by flux ratio anomalies (also see Evans & 
Witt 2002; Quadri et al. 2003; Moller et al. 2003). Plau- 
sibility arguments and astrophysical expectations might 
be invoked to favor one possibility over another. For ex- 
ample, stars are known and dark matter clumps expected 
to be abundant in lens galaxies, while high-order multi- 
pole modes are not, so one might prefer to attribute flux 
ratio anomalies to stars and dark matter clumps. But it 
is important to acknowledge the prejudices inherent in 
such arguments. 

Fortunately, there are excellent prospects for obtaining 
data that move beyond broad-band fluxes to determine 
the nature of the small-scale structure. One possibil- 
ity is to look for time variability that is an unmistak- 
able signature of microlcnsing by stars (e.g., Wozniak et 
al. 2000; Schechter et al. 2003). A second good possi- 
bility is to look for Einstein ring images of the quasar 
host galaxy in deep near-IR images. The Einstein rings 
would vastly improve the constraints on the global prop- 
erties of the lens potential and reveal whether high-order 
multipole modes are acceptable (Kochanek et al. 2001; 
Kochanek & Dalai 2003). A third possibility is to use an 
aspect of the cusp relation that we have ignored, namely 
the sign of R CU s P - Several recent studies have indicated 
that localized structure is sensitive to the image parities, 
most often suppressing negative-parity images and oc- 
casionally increasing the magnification of positive-parity 
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images (Schechter & Wambsganss 2002; Keeton 2003), 
while global modes make no distinction between differ- 
ent images. Thus, in an ensemble of lenses with flux ra- 
tio anomalies, the presence or absence of skewness in the 
set of R cusp values could reveal whether the small-scale 
structure is local or global. Indeed, Kochanek & Dalai 
(2003) already find evidence that many flux ratio anoma- 
lies can be attributed to suppression of negative-parity 
images, which suggests that many flux ratio anomalies 
are in fact caused by milli-lensing and/or micro-lensing. 
Finally, resolved high-resolution spectroscopy of systems 
with flux ratio anomalies offers the intriguing possibility 
of determining the physical scale of the structure (Mous- 
takas & Metcalf 2003; Wisotzki et al. 2003). Thus, the 
future is very bright for using detailed study of lenses 
with flux ratio anomalies to learn about small-scale struc- 
ture in lens galaxies. Since all of these applications begin 
with the identification of flux ratio anomalies, the cusp 
relation analysis will be of fundamental importance to all 
of them. 
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APPENDIX 

A. UNIVERSAL RELATIONS FOR CUSPS 

In this Appendix we study the general properties of the lensing map near a cusp catastrophe to derive generic 
relations between the image positions and magnifications that should be satisfied whenever the source is sufficiently 
close to a cusp. This analysis applies to ordinary cusps, in which the two branches of the curve approach the cusp from 
opposite sides of the line that is tangent to the cusp. The analysis may not be valid for ramphoid cusps, an alternate 
situation where the two branches approach the cusp from the same side of the tangent line (for examples, see Petters 
& Wicklin 1995; Oguri et al. 2003). Only ordinary cusps have been observed, and ramphoid cusps are expected to be 
rare in lensing situations of astrophysical interest. 

A.l. Local orthogonal coordinates 

Consider the lens equation 

y = x — grad ip (x) . (Al) 

Assume that the induced lensing map, rj(x) = x — grad?/>(x), from the lens plane L to the light source plane S is 
locally stable, which yields that the caustics of r\ are either folds or cusps (Petters et al. 2001, p. 294). Then translate 
coordinates in the lens and light source planes so that the cusp point of interest is at the origin and the critical point 
mapping to the cusp is also at the origin, i.e., rj(0) = 0. By abuse of notation, the resulting translations of x and y 
will still be denoted by those symbols. 

We now define a change of coordinates about the origins of the lens and light source planes: x — * 8, y — > u. The 
coefficients of the quadratic terms of the Taylor expansion of ip at the origin are 

a=^'n(0), S = Vi2(0), c=^ 22 (0), (A2) 

where the subscripts indicate partial derivatives relative to x = (xi,x%). Since the origin is a cusp critical point, 
(1 — 2d) and (1 — 2c) cannot both vanish (Petters et al. 2001, p. 349). Without loss of generality, we shall assume that 
(1 — 2d) 7^ 0. Define an orthogonal matrix as follows: 



M = 



1 



1 - 2d -b 
b 1 -2a 



(A3) 



'(l-2d) 2 +6 2 

The new coordinate systems arc defined by 

e = {0 u e 2 ) = Mx, u = («i,«a) = My. (A4) 

Note that this is the same coordinate change in the lens and light source planes, and the coordinates depend on the 
potential. 

It can be proven rigorously that the lensing map r\ in the orthogonal coordinates (|A4|) can be approximated in a 
neighborhood of the cusp critical point at the origin by a simple polynomial mapping (Petters et al. 2001, pp. 341-353; 
also see Schneider et al. 1992, p. 193): 

ui = c 0i + - 6l , u 2 = b9 1 Q 2 +ae%, (A5) 
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where 

a = --^2222(0) , b — -^122(0) 7^ , 

c = 1-Vh (0)^0, 2ac-6 2 ^0. (A6) 
The partial derivatives of ip are relative to the original coordinates x = (xx,X2)- The origin in the light source plane 
is called a positive cusp if 2ac > b 2 and a negative cusp if lac < b 2 . A source inside a positive cusp has, locally, two 
images with positive parity and one with negative parity; the reverse is true for negative cusps. 

A. 2. Position relations 

Using the lens equation l|A5J) . the three local lensed images associated with a source inside and close to the cusp 
have the following positions (e.g., Gaudi & Petters 2002): 

e *=Cj--b*i > **)> i= 1,2,3. (A7) 

The Zi are the three real solutions of the cubic equation 

z 3 + pz + q = , (A8) 

where 

2b 2c . 

P= 2^c—V 2Ul=PU1 ' q= ~2a^ U2 = - qU2 - (A9) 
Note that when the source is inside the cusp the discriminant, 

is negative so eq. I|A8(I does have three real roots. 
The usual factoring of a cubic polynomial yields: 

= (z-z 1 ){z-z 2 )(z-z 3 ), (All) 

= z 3 - [z\ + z 2 + z 3 ] z 2 + \z\Zi + ziz 3 + z 2 z 3 \ z ~ [ziz 2 z 3 ] . (A12) 
Identifying coefficients with eq. I|A8(1 yields three relations between the image positions: 

zi + z 2 + z 3 = 0, (A13) 

ziz 2 + ziz 3 + z 2 Z3=pui, (A14) 

ziz 2 z 3 = qu 2 . (A15) 
These are universal relations satisfied by the image positions of a triplet associated with a source near a cusp. Two 

addit ional relations can be obtained respectively by squaring HA13fl and using (|A14(i . and squaring l|A14|l and using 

EE): 

z\ + z\ + z| = -2pui, (A16) 
( Zl z 2 ) 2 + ( Zl z 3 ) 2 + (z 2 z 3 ) 2 = (-p Ul ) 2 . (A17) 
These relations are not independent of ilA13|) - (|A15|) . but they are more useful in certain circumstances (as seen below). 

A. 3. Magnification relations 
The signed magnification of each image Oi in the triplet associated with the cusp is given by 

_1 = P_ 

det[Jacu](0;) b(pui + 3 z 2 
where Jacu is the Jacobian matrix of the lensing map l)A5|l . Note that Jacw = M Jacy, so with M an orthogonal 
matrix we verify that the magnification is independent of our choice of coordinates: det [ J ac it] = det[Jacy]. 

Three known universal relations between the magnifications \ii are as follows (Schneider & Weiss 1992; Zakharov 
1995; Petters et al. 2001, p. 339): 

Mi + M2 + M3 = 0, (A19) 
M1M2 + M1M3 + A^3 = - 36b2£) , (A20) 

™ 3= iol^' (A21) 

where D is giv e n by eq. IjAlOf) . These relations can be verified by direct calculation from l|A18|) . using the position 
relations (|A15|I — (|A1T|) for simplifications. In analogy with the position relations, we can derive additional magnification 
relations: 

A 3 

2 1 2 1 2 P M l /a oT\ 

Mi + M2 + th = Tg-^FTJ ' ( A22 ) 

(M1M2) 2 + (MiM3) 2 + (^M3) 2 =(^^) • (A23) 
These quadratic magnification sum rules have not appeared in the literature before. 



V* = ^n^..ua ^ = U^T. o 2, . * = 1.2,3, (A18) 
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B. SIMPLE LENS POTENTIALS 



In this Appendix we derive exact solutions to the lens equation to use as a benchmark for understanding the cusp 
relations. Exact solutions are possible only for certain lens potentials, and then only for sources on a symmetry axis. 
We consider two families of potentials: a spherical galaxy with a power law density profile plus an external shear; and 
a singular isothermal ellipsoid with multipole density perturbations plus an external shear aligned with the major or 
minor axis of the galaxy. 



Consider the lens potential 



B.l. Power law galaxy with shear 



^(r,^) = -R 2 c r a r a - Jr 2 cos20. 
a 2 



The first term represents a spherical galaxy with a power law profile for the surface mass density, 

2-a 

K(r) 



S(r) 



-'crit 



Re 



(Bl) 



(B2) 



where -Rein is the Einstein radius. The case a = 1 corresponds to a singular isothermal sphere (SIS), while the 
cases a < 1 and a > 1 correspond respectively to steeper and shallower profiles, respectively. The second term in the 
potential represents an external tidal shear with amplitude 7. Without loss of generality, we are working in coordinates 
such that the shear is aligned with the horizontal axis (7 > 0) or the vertical axis (7 < 0). 

Using polar coordinates in the image plane and Cartesian coordinates in the source plane, the lens equation has the 
form 



Vi- 



- r cos < 



1+7- 



2-a 



yi = r sin < 



R c 



and the lensing magnification /1 is given by 



1 - 7 2 - (1 - a) 



Rc\\ 



4-2a 



2-a 



[a + (2 — a)"f cos 2(f>] . 



(B3) 
(B4) 

(B5) 



The critical curve in the image plane is the curve where = 0, and it maps to the caustic in the source plane. The 
caustic has a cusp on the horizontal axis at position (yi c ,0), which corresponds to a point on the critical curve at 
position (a;i c ,0), where 

Rcin (B6) 



X\c = 



!Jlc 



(1 - 7 )V(2-a) ' 
27-Rein 



(B7) 



(1 - 7)V(2-«) ' 

The Taylor series coefficients used to define the local orthogonal coordinate system in Appendix IA A. H are as follows 

1 

C= 2 ' 



: i [_l + a (l-7)] 



0. 



2-a 

'ein 



2R 
2 



(1-7) 



(1-7) 



(4-a)/(2-a) 



(3-a)/(2-a) 



c = 2 



P = 



a(l - 7) . 



7(1 - 7)V( 2 - Q ) 



(B8) 
(B9) 

(BIO) 



Hence the transformation matrix M in eq. ljA4fl is the identity matrix, so the 9 and u coordinate systems are simply 
the x and y coordinate systems translated so the cusp point is at the origin. 

Note that although the potential %p is not well defined in the limit a — > 0, the lens equation and magnification and 
other quantities are perfectly well defined and correspond to a point mass in a shear field. Furthermore, in this limit 
the surface mass density S is a (5-function as expected for a point mass. Hence in this formalism we can consider the 
case a = to correspond to a point mass lens. 

Consider a source on the horizontal axis inside the caustic; for 7 > (7 < 0) this correspond to the major (minor) 
axis of the lens potential. The lens equation can be solved exactly because of symmetry. There is at least one image on 
the xi-axis, 14 and two images off the xi-axis. By symmetry, the two off-axis images are identical modulo some signs. 

14 There may or may not be an image on the xi-axis on the opposite side of the origin from the source, depending on whether the cusp 
is "clothed" or "naked" (e.g., Schneider et al. 1992; Petters et al. 2001). We are interested only in the image on the xi-axis on the same 
side of the origin as the source. 



17 



To find the positions of these two images, note that with j/2 = and <p ^ the only way for eq. (|B4|) to be satisfied is 
for the t erm in square brackets to vanish. This condition yields the polar radius, which can then be substituted into 
cq. (|B3ll to find the polar angle. Thus, the positions of the two off-axis images, which we label A and C, are 



R , 



rA = r c 



(1 - 7 ) 1 /(2-a) 



COS 



Vic 



where y\ c is given by eq. i|B7|) . Their magnifications of these two images are 



^=Mc = ^2 7 (l- 7 )(2-a) 



1-1 J* 

Vic 



(Bll) 
(B12) 

(B13) 



The image separation for this triplet is simply d = 2r J 4sin0 J 4. 

For the on-axis image, which we label B, eq. I|B4|) is satisfied trivially (y2 = and <\>b = 0). Eq. I|B3|) can be solved 
analytically for integer and half-integer values of a, yielding: 



= 

1 
2 



a= 2 



rB = 



TB 



yi + ^y'i + 4(l + j)R 



2 

ein 



2(1+7) 



(see below) 



(e + yi) 2 

3^(1 + 7) 

Vl + -Rein 

1+7 

i? ei „ + 2(1 + 7)2/1 + y/R~[R c 



4(l+7)i/i 



In the result for a 



1/2, £ is given by 

e 3 = f (i+7)4 



vl 



2(1 + 7) 2 



{3(1 + 7 )i? e 3 in [27(1 + 7 )i? c 3 in - %?] } 



1/2 



(B14) 
(B15) 
(B16) 
(B17) 

(B18) 



The magnification /is of image B can then be computed from eq. l|B5ll . 

This analysis applies only to sources on the symmetry axis of the lens, but we can begin to understand what happens 
when the source is moved off-axis by examining derivatives with respect to yi- The first derivative of i? C usp vanishes 
by symmetry, 



OR, 



cusp 



dy2 



= 0. 



(B19) 



2/2=0 



so the axis is a local extremum. The second derivative, which determines whether it is a local maximum or minimum, 
can be computed explicitly for an SIS plus shear potential. After lengthy but straightforward algebra, we find 



d 2 R 



cusp 



dy. 



(1+7) 2 



2/2=0 



? • 4 

7 Z sin 



(6»/2) [3 + 2 7 + (1 + 7) cos(0/2) + 7 cos 9} 2 
4 7 (1 - 7) + (7 + 67 - 7 2 ) cos(0/2) + 167(2 + 7) cos 2 (0/2) 

+ (5 + I87 + 13 7 2 ) cos 3 (0/2) + 12 7 (1 + 3 7 ) cos 4 ((9/2) 



(B20) 



The factor on the first line is manifestly negative, while the quantity in square brackets on the second and third lines 
is positive over the entire interesting range < 8 < it and | 7 | < 1. Thus, the second derivative is negative, and hence 
■Rcusp is a maximum on the axis. While this proof formally holds only for the SIS plus shear potential, intuition and 
Monte Carlo simulations suggest that it is not restricted to this model. On-axis sources therefore provide a simple and 
important upper bound on i? C usp- 

B.2. Generalized "isothermal" galaxy with shear 
Consider the potential/density pair 



if)(r, ft) = rF(<p) - ^r 2 cos 2(f> . 



k(V, (j>) - 



G(4>) 
2r 



(B21) 
(B22) 
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where, from the Poisson equation, F(<fi) and G(<fi) are related by 

G(<f>)=F(<f>)+F"(cl>). (B23) 

The density and the hrst term in the potential correspond to a mass distribution that is scale-free in the radial direction 
and produces a flat rotation curve; such a model is often referred to as "isothermal" in the lensing literature. The mass 
distribution is allowed to have an arbitrary angular shape specified by the functions F(cj)) and G(</>). This family of 
models includes both the singular isothermal ellipsoid and the singular isothermal elliptical potential but is much more 
general, and its lensing properties have been studied by Witt, Mao, & Keeton (2000), Evans & Witt (2001, 2002), and 
Zhao & Pronk (2001). The second term represents an external tidal shear with amplitude 7, in coordinates such that 
the shear is aligned with the horizontal axis (7 > 0) or the vertical axis (7 < 0). 

In order to make analytic progress with this model, we assume that the shape function F((f>) is an even function, 
i.e., F((f>) = F(~<fi). In other words, we assume that the galaxy is symmetric about the horizontal axis. The shear 
we consider therefore does not have an arbitrary orientation, but is either aligned with or orthogonal to the galaxy's 
symmetry axis. Although not completely general, these two cases should bound the interesting range of shears. 

Using polar coordinates in the image plane and Cartesian coordinates in the source plane, the lens equation has the 
form 

y x = (1 + 7) r cos <j) - F{4>) cos cf> + F'{<j>) sin , (B24) 
Vi = (1 - 7) r sin cf> - F{4>) sin - F'(<j>) cos (f> , (B25) 
and the lensing magnification is \i given by 15 

/i" 1 = 1 - 7 2 - 2(1 + 7 cos 2<j>) «(r, cj>) . (B26) 
The critical curve can be written in parametric form as 

1 + 7 cos 2(f> 

which can be used in the lens equation to obtain a parametric expression for the caustic. The condition F(<fi) = F(—<fi) 
ensures that there is always a cusp at = 0, whose location in the image and source planes is 

x lc = ^, (B28) 
1-7 

2 7 F(0) + (1 + 7)F"(0) 
Vlc = : • (B29) 

1-7 

In general this cusp is a simple cusp, but for some combinations of the shape function G{0) and shear 7 it can be 
part of a higher-order butterfly catastrophe. We find that butterfly catastrophes are rare on the major axis of the 
lens potential, but can be relatively common on the minor axis when the potential has significant power in high-order 
multipole modes. 

The Taylor series coefficients used to define the local orthogonal coordinate system in Appendix IA A. H are as follows: 

£ -0, c=-, (B30) 



rc(fl= Y (B27) 



2 ' ' 1 

(I-7) 3 3G(0)-G"(0) 



6 G(0) 



3 



c=l + 7 , (B31) 



h= (1-7) 2 fi= 6G^! fmol 

G(0) ' P (1_ 7 )[6 7 G(0)-(1+ 7 )G"(0)] ' 1 ] 

Hence the transformation matrix M in eq. ljA4fl is the identity matrix, so the 6 and u coordinate systems are simply 
the x and y coordinate systems translated so the cusp point is at the origin. 

A source on the horizontal axis inside the caustic has at least one image on the xi-axis and two images off the 
si-axis. Because of the reflection symmetry F{4>) = F (—(/>), the two off-axis images are id entical modulo some signs. 
The polar radius for these two images, labeled A and C, is found by requiring that eq. (|B25(I have a non-trivial solution 
(i.e., 0^0): 

r A = r c = - J— [F(<t> A ) + F'((t> A ) cot <j> A \ . (B33) 
1-7 

Their polar angles satisfy 4> A — —<frc — 0/2 where 8 is the openin g ang le defined in <J5] The image separation for this 
triplet is simply d — 2r A sin cf> A . The source position is, from eq. (|B24|) . 

Vi = -7 ]. 7 hF{0 A )sm2<f> A + (l+ 1 cos2<f> A )F'(<f> A )}. (B34) 

(1- 7) sin 0,4 

15 Note that in the absence of shear (7 = 0), the magnification is simply fj, = (1 — 2k)~ 1 and the critical curve is the isodensity contour 
K = 1/2 of the galaxy (Witt et al. 2000; Evans & Witt 2001, 2002; Zhao & Pronk 2001). 



19 




Finally, the position of the on-axis image (labeled B) is found by solving eq. (|B24|) with ips = 0: 

rs - (B35, 



Using eq. (|B26(I . we find the magnifications of the three images to be 

a? - - 1 - t - a - F( ^ ) ;,I',t^, • < b ») 

-1 n 2 ,, , , 2 F(0)+F"(0) 
M fl =1-7 -(1 + 7) F(0)+yi ■ (B37) 

A specific case of interest is a singular isothermal ellipsoid (SIE), which has the shape functions (Kassiola & Kovner 
1993; Kormann, Schneider, & Bartclmann 1994; Keeton & Kochanek 1998) 

G si M= n = m , (B38) 
V 1 — £ cos 20 



F sic (0) = /,>( 



2e cos0 \ . , _i / V2e sin0 
-^^^=^^^= + sin tann — ^^=^^^= 
Vl — £ cos 20 / \ VI — £ cos 20 



(B39) 



(The parameter e is related to the minor-to-major axis ratio g of the ellipse by £ = (1 — q 2 )/(l + q 2 ), and it is a 
convenient parameter for this formalism; however, in the main body of the paper we always quote the true ellipticity 
e = 1 — q.) When e > (e < 0) this formalism describes a source on the major (minor) axis of the galaxy. When 
£ and 7 have identical (opposite) signs, the shear is aligned with (orthogonal to) the galaxy's major axis. We note 
that when thinking in terms of a multipole expansion, the SIE has power in all even multipole moments; the multipole 
coefficients, 

i? cin r cos wo d ^ 



27r J y/l — s cos 20 
are shown in Figure IBlOl 

We also consider adding perturbations that represent departures from elliptical symmetry in the density distribution. 
Examples of such perturbations are "boxy" or "disky" isophotes, or disk-like components, all of which are observed 
(Bender et al. 1989; Saglia et al. 1993; Kelson et al. 2000; Rest et al. 2001; Tran et al. 2003) and predicted (Heyl et al. 
1994; Naab & Burkert 2003; Burkert & Naab 2003) in early-type galaxies. It is convenient to express the perturbations 
in terms of multipole modes. The shape functions for an m-th order mode are 

G m (0) = aP l ert cos(m0), (B41) 

port 

F^)^ cos(m0), (B42) 
1 — m 

where the perturbation amplitude a„ rt ^ s defined such that the deviation of an isodensity contour (say, k = 1/2 
although since the potential is scale-free the choice is irrelevant) from a pure ellipse is 

Sr = a P° rt cos(m0) . (B43) 
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This definition is exactly equivalent to the amplitude used to quantify isophotes shapes in observed galaxies (e.g., 
Bender et al. 1989; Saglia et al. 1993; Rest et al. 2001). It is closely related to the fractional change in the surface 
density, 

X K „pcrt 

— = cos(m0) VT^I^oi^ . (B44) 

The perturbation amplitude has dimensions of length, but we follow observational studies and normalize it by the size 
of the reference ellipse. We can therefore adopt <24 perturbation amplitudes directly from the observational studies. 
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Table Bl. Observed 4-Image Lenses 



Lens and i?cusp 



References 


Tvne 


/"?,-,;.-. ( i 


Triplet 


8 (°) 


d (") 


Data 


Model 


B0128+437 


radio, 


0.20 


BCD 


236.7 


0.50 


0.30+0.06 


0.40 


(13) 


fold 




ACD 


197.4 


0.55 


0.52+0.05 


0.44 






ABD 


123.3 


0.34 


0.01+0.06 


0.03 








ABC 


162.6 


0.55 


0.51+0.05 


0.49 


HE 0230-2130 


optical, 


0.83 


A 2 BC 


197.2 


1.66 


0.35+0.06 


0.09 


(1) 


fold 




AiBC 


231.6 


2.19 


0.84+0.02 


0.77 






AiA 2 C 


162.8 


1.66 


0.03+0.07 


0.01 








AiA 2 B 


128.4 


2.19 


0.28+0.06 


0.41 


MG 0414+0534 


near-IR, 


1.08 


A 2 BC 


216.0 


2.13 


0.40+0.05 


0.57 


(1) 


fold 




AiBC 


258.5 


2.13 


0.77+0.03 


0.78 






A1A2C 


144.0 


2.08 


0.04+0.06 


0.09 








A1A2B 


101.5 


2.03 


0.31+0.06 


0.12 


HE 0435-1223 


optical, 


1.18 


BCD 


179.0 


2.25 


0.34+0.05 


0.30 


(16) 


cross 




ACD 


201.7 


2.56 


0.48+0.05 


0.50 






ABD 


181.0 


2.25 


0.08+0.06 


0.20 








ABC 


158.3 


2.56 


0.43+0.05 


0.34 


B0712+472 


radio, 


0.65 


BCD 


200.4 


1.27 


0.33+0.06 


0.62 


(8, 9) 


cusp 




ACD 


283.1 


1.25 


0.89+0.01 


0.91 






ABD 


159.6 


1.27 


0.06+0.07 


0.07 








ABC* 


76.9 


1.05 


0.26+0.06 


0.08 


B0712+472 


opt/IR, 


0.65 


BCD 


200.4 


1.27 


0.44+0.05 


0.62 


(1,9) 


cusp 




ACD 


283.1 


1.25 


0.69+0.03 


0.91 






ABD 


159.6 


1.27 


0.32+0.06 


0.07 








ABC* 


76.9 


1.05 


0.60+0.04 


0.08 


RX J0911+0551 


near-IR, 


0.96 


A2A3B 


179.4 


3.26 


0.45+0.05 


0.44 


(1) 


cusp 




A1A3B 


290.4 


3.08 


0.59+0.04 


0.67 








AiA 2 B 


180.6 


3.26 


0.14+0.06 


0.38 








A1A2A3* 


69.6 


0.96 


0.23+0.06 


0.00 


SDSS J0924+0219 


optical, 


0.87 


BCD 


217.7 


1.61 


0.04+0.08 


0.46 


(7) 


cross 




ACD 


142.3 


1.61 


0.36+0.07 


0.09 






ABD 


156.4 


1.79 


0.89+0.02 


0.29 








ABC 


203.6 


1.79 


0.58+0.05 


0.56 


PG 1115+080 


optical, 


1.14 


A 2 BC 


233.3 


2.16 


0.52+0.05 


0.59 


(6) 


fold 




AiBC 


218.8 


2.43 


0.76+0.03 


0.70 






AiA 2 C 


141.2 


2.43 


0.32+0.06 


0.18 








A1A2B 


126.7 


1.86 


0.10+0.06 


0.06 


1RXS J1131-1231 


optical, 


1.81 


BCD 


290.5 


3.18 


0.87+0.01 


0.90 


(15) 


cusp 




ACD 


181.0 


3.20 


0.34+0.06 


0.31 






ABD 


179.0 


3.20 


0.15+0.06 


0.29 








ABC* 


69.5 


2.38 


0.35+0.06 


0.06 


HST 12531-2914 


optical, 


0.54 


BCD 


172.4 


1.04 


0.39+0.05 


0.28 


(1, 14) 


cross 




ACD 


187.6 


1.04 


0.46+0.05 


0.14 






ABD 


206.8 


1.36 


0.31+0.06 


0.59 








ABC 


153.2 


1.36 


0.17+0.06 


0.33 
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Table Bl. Observed 4-Image Lenses — Continued 



Lens and 












ficusp 




References 


Type 


Rein (") 


Triplet 


e (°) 


d(") 


Data 


Model 




HST 14113+5211 


optical, 


0.83 


BCD 


168.9 


1.35 


0.14+0.06 


0.09 


(4) 


cross 




ACD 


161.3 


2.28 


0.36+0.05 


0.42 






ABD 


191.1 


1.42 


0.15+0.06 


0.16 








ABC 


198.7 


2.28 


0.67+0.03 


0.63 


H1413+117 


near-IR, 


0.56 


BCD 


198.6 


1.35 


0.48+0.05 


0.28 


(1) 


cross 




ACD 


186.2 


1.10 


0.37+0.05 


0.83 






ABD 


173.8 


1.10 


0.26+0.06 


0.78 








ABC 


161.4 


1.35 


0.24+0.06 


0.61 


HST 14176+5226 


optical, 


1.33 


BCD 


172.8 


2.36 


0.04+0.06 


0.03 


(1. 14) 


cross 




ACD 


198.1 


3.26 


0.54+0.04 


0.64 






ABD 


187.2 


2.36 


0.13+0.06 


0.08 








ABC 


161.9 


3.26 


0.54+0.04 


0.50 


B1422+231 


radio, 


0.76 


BCD 


187.2 


1.29 


0.35+0.06 


0.35 


(5, 12) 


cusp 




ACD 


283.0 


1.29 


0.96+0.01 


0.94 








ABD 


172.8 


1.25 


0.05+0.07 


0.15 








ABC* 


77.0 


1.29 


0.18+0.06 


0.12 


B1555+375 


radio, 


0.23 


BCD 


209.3 


0.42 


0.14+0.07 


0.56 


(11) 


fold 




ACD 


257.4 


0.42 


0.90+0.01 


0.89 






ABD 


150.7 


0.42 


0.21+0.06 


0.03 








ABC 


102.6 


0.41 


0.45+0.05 


0.14 


B1608+656 


radio, 


0.72 


BCD 


191.5 


2.04 


0.16+0.06 


0.06 


(10) 


fold 




ACD 


168.5 


2.04 


0.19+0.06 


0.24 






ABD 


261.0 


2.10 


0.79+0.02 


0.89 








ABC 


99.0 


2.10 


0.49+0.05 


0.49 


B1933+503 


radio, 


0.49 


3,4,6 


143.0 


0.82 


0.39+0.05 


0.01 


(2) 


fold 




1,4,6 


199.7 


1.16 


0.70+0.03 


0.63 






1,3,6 


217.0 


0.91 


0.21+0.06 


0.29 








1,3,4 


160.3 


1.16 


0.72+0.03 


0.42 


B2045+265 


radio, 


1.13 


BCD 


183.9 


1.93 


0.05+0.06 


0.49 


(3) 


cusp 




ACD 


325.1 


1.93 


0.88+0.01 


0.98 






ABD 


176.1 


1.92 


0.21+0.06 


0.19 








ACD* 


34.9 


0.84 


0.52+0.04 


0.02 


Q2237+030 


optical, 


0.85 


BCD 


186.5 


1.65 


0.29+0.06 


0.18 


(1) 


cross 




ACD 


173.5 


1.65 


0.20+0.06 


0.12 






ABD 


146.2 


1.83 


0.71+0.03 


0.39 








ABC 


213.8 


1.83 


0.52+0.05 


0.62 



Note. — Results for image triplets in the nineteen published 4-image lenses. Column 2 gives the image configuration (fold, cusp, or 
cross) and indicates whether the flux ratios are measured at optical, near-IR, or radio wavelengths. The uncertainties in the observed values 
of -Rcusp are obtained by assuming 10% uncertainties in the image fluxes; see JH]for more discussion. The predicted values of i? C usp are 
computed with standard lens models. For the cusp lenses B0712+472, RX J0911+0551, 1RXS J1131-1231, B1422+231, and B2045+265, 
the cusp image triplet is indicated by *. Note that B0712+472 appears twice because we report data from both radio and optical/near-IR 
wavelengths. The references are as follows: (1) CASTLES (see http://cfa-www.harvard.edu/castles); (2) Cohn et al. 2001; (3) Fassnacht 
et al. 1999; (4) Fischer et al. 1998; (5) Impey et al. 1996; (6) Impey et al. 1998; (7) Inada et al. 2003; (8) Jackson et al. 1998; (9) Jackson 
et al. 2000; (10) Koopmans & Fassnacht 1999; (11) Marlow et al. 1999; (12) Patnaik et al. 1999; (13) Phillips et al. 2000; (14) Ratnatunga 
et al. 1995; (15) Sluse et al. 2003; (16) Wistoski et al. 2002. 



